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NOMENCLATURE 


eddy viscosity constants 

term in Van Driest damping factor 

2 

area of the injection device, ft 

density-viscosity product ratio, py/(p„y_) 

mass fraction of the freestream species 

skin-friction coefficient based on freestream conditions 

mass fraction of species i 

specific heat of perfect gas, mV(sec -°K) 

2 

binary diffusion coefficient, m /sec 

tangential velocity ratio, u/U^ 

form factor in multi-layer eddy viscosity model 

stagnation enthalpy ratio, H/H^ 

2 2 

Stagnation enthalpy, m /sec 

2 2 

Static enthalpy, m /sec 
eddy viscosity constant, 0.4 
reference length, m or cm 
molecular Lewis number 
turbulent Lewis number 
mixing length, m or cm 
freestream mixing length 
mixing region mixing length 
slot mixing length 
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Mach number 

2 

pressure, N/m 

molecular Prandtl number 

turbulent Prandtl number 

gas constant 

Reynolds number 

slot height, m or cm 

temperature, °K 

slot lip thickness, m or cm 

tangential velocity component, m/sec 

friction velocity, 

tangential velocity in law of the wall variables, u/u^ 

transformed normal velocity 

normal velocity component, m/sec 

width of the mixing region 

weighting factor 

distance along surface, m or cm 

distance normal to surface, m or cm 

normal coordinate in law of the wall variables, y u^/v 

y at Z = 0.99 

y at Z = 0.5 

y at Z = 0.01 

mass fraction, C^/C^. 
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Normalized mass fraction ratio, (C.. - C. .,)/C. _ - C. 

• ifC I »w 

Or dUo 

pressure gradient parameter, ^ ~g|’ 

Klebanoff i ntermi ttency factor 
boundary-layer thickness 
boundary-layer thickness on flat plate 
6o + t + s 

boundary-layer displacement thickness 
incompressible boundary-layer displacement thickness 
eddy viscosity 

e/y 

* * * * - 1/2 
Van Dyke parameter (p^ L 

Levy-Lees normal coordinate 

incompressible boundary- layer momentum thickness 
viscosity 

Levy-Lees streamwise coordinate 
density 

shear stress at wall, v ^ 

L J 

Subscripts and Superscripts 

condition at outer edge of boundary layer 
species i 

★ 

denotes limits of integration for and from y^ to 6 
condition at wall boundary 


stagnation condition, flat plate conditions 
differentiation with respect to 5 
freestream condition 
differentiation with respect to n 
dimensional quantity 


injectant conditions 



I. INTRODUCTION 


In recent years, injection of a fluid into the boundary layer has 
been studied both numerically and experimentally as a technique for ther- 
mal protection and skin friction reduction. Two of the most common injec- 
tion configurations studied have been the porous wall and the tangential 
slot. In the former configuration, air or a foreign gas is injected into 
the boundary layer through a porous section of the wall. The tangential 
slot consists of a step in the wall through which the injectant enters 
the boundary layer parallel to the frees tream flow. 

Spalding and Patankar (Ref. 1) were able to obtain a solution for 
the tangential slot by using two regions for the flow field. In the 
first region, the flow was treated as a mixing of two planar flows and 
the wall region of the flow was neglected. When the mixing region had 
spread to the wall, the standard solution procedure for a turbulent bound- 
ary layer was used. With this model, Spalding and Patankar obtained good 
predictions of the effects of slot injection at some distance downstream, 
but predictions of the effects in regions near the slot were not as good. 

Beckwith and Bushnell (Ref. 2) developed a finite-difference method 
which could be used near the slot as well as far downstream. This method 
used a modified form of the eddy viscosity expression developed by Bushnell 
(Ref. 3). In this expression, the Prandtl mixing length in the region near 
the slot was adjusted to account for the effects of slot injection. The 
mixing length in the mixing region downstream of the slot lip was propor- 
tioned to the width of the mixing region and thus increased downstream 
as the region spread between the slot and freestream flows. By adding 
the species conservation equation to the system of governing equations 
(Ref. 2), the growth of the mixing region was calculated and used to com- 
pute the mixing length. The method of Ref. 2 was limited to the case where 
the injected species was the same as the external flow. 


Miner and Lewis (Ref. 4) presented a finite-difference method for 
predicting compressible, turbulent boundary-layer flows with tangential 
slot injection. Inclusion of the species conservation equation allow- 
ed predictions of the mixing between the slot and freestream flows and 
use of the Beckwith-Bushnell eddy viscosity model. Comparisons with the 
experimental data of Kenworthy and Schetz (Ref. 5) were made for three 
eddy viscosity models to test the accuracy of the predictions in a re- 
gion near the slot. The method also contained a procedure to predict 
the effects of the interaction between the boundary-layer displacement 
thickness and the external pressure field. 

In a later work (Ref. 6), Miner and Lewis extended their method to 
include non-reacting foreign gas injection. With a foreign gas as the 
injectant, the single speciesgas property model had to be replaced by a 
binary gas model. The properties used in the binary gas property model 
are given in Ref. 6. 

Experimental studies of this problem have been conducted by a number 
of authors. McRee,et al . (Ref. 7) considered the effects of normal injec- 
tion and slot injection on the skin friction at Mach 3. Cary and Hefner 
(Refs. 8, 9) studied the effects of tangential -slot injection on the film- 
cooling effectiveness and skin friction in a hypersonic flow of Mach 6. 
Finally, Kenworthy and Schetz (Ref. 5) conducted an experimental study of 
slot injection into a Mach 2.4 freestream. 

The present study used a finite difference method to compare the 
effectiveness of slot injection, porous injection, and a simple combina- 
tion of the two in reducing the skin friction. Comparisons are made with 

experimental data from a similar study by Schetz and Van Overeem (Ref. 10). 

2 

The conditions for the experimental study were Mach 2.9, = 69.9 N/cm 

(100 psia), and T^ * 294 °K (530 °R). The data taken for comparison in- 
cluded wall -pressure distributions, Mach number profiles and pitot pres- 
sure profiles at four axial stations and wall shear as measured by a 
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floating element balance. Wall shear data were also taken by two Preston 
tubes with diameters of 0.241 cm and 0.073 cm. 

The finite-difference method developed by Miner and Lewis (Ref. 4) 
was used in the present work. Calculations of the skin friction with and 
without pressure interaction effects were made using the two-layer eddy 
viscosity model. For the slot and slot/porous configurations, the multi- 
layer eddy viscosity model of Beckwith and Bushnell was also used for com- 
parisons. Mach number profile calculations downstream of the slot were 
made and compared with the experimental data. 

Finally, the finite-difference method was used to study the effect 
of combined slot-porous injection configurations on reducing the skin 
friction. Effects of pressure interaction and the differences in the two 
eddy viscosity laws were also considered. 
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II. ANALYSIS 


The equations for two-dimensional, turbulent, boundary-layer flows 
of a non-reacting, two-component mixture of gases are presented in non- 
dimensional, Levy-Lees variables. These equations are then presented 
in the parabolic form for solution by the finite-difference method of 
Anderson and Lewis (Ref. 11). Specification of the fluid properties, 
eddy viscosity model, and a brief summary of the global pressure itera- 
tion scheme is included. 


2.1 Governing Equations 

The boundary layer equations were first non-dimensional ized as pro- 
posed by Van Dyke and then transformed using Levy-Lees variables. A ref- 
erence temperature and reference viscosity were defined as 

<’) 

Vef - (<ef) (2) 

★ 

Using any suitable reference length L the Van Dyke parameter became 

9 * * * ★ 

'VD ' Vef''=X'- <3) 

The non-dimensional variables which were used in the present method 
thus become 

X = x*/L* (4) 

it ^ 

y ^ y ^vD 

P * P /pX (6) 
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* * 

p * p 


★ ★ 
T = T /T 


ref 


* * 
U * U /U 


* * 

V * V /U E 


VD 


(7) 

(8) 

(9) 

( 10 ) 


The resulting boundary-layer equations were transformed using the Levy- 
Lees variables n and c defined by 


dS = p_Uq)j_ dx 
e e e 


( 11 ) 


and 


pUp 

dn = -7==:dy 


ViF 

The turbulent boundary layer equations were then expressed as: 
Continuity : 

2sF^+ V + F = 0 


Momentum: 


Energy : 

2 cFg^ + Vg' = 


25FF^ + VF' = e(pg/p-F^) + [c(l + e"") F' ] 


^ (1 e^)g' 


Pr 
Pr. 


C 

Pr 


(Le . Le, e^) 


( 12 ) 


(13) 


(14) 


C 

Pr 


(15) 
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Species : 

2cFZj + Vr .[[^(Ut^Lej e^)l 2' 1 ' 

I L t J / 

where V Is the transformed normal velocity component 

2L 


V = 


Q U u 

e^e 


^ , (py tgQ 

L " 


(16) 


(17) 


and where p'v' is the time average of the product of the fluctuating 
density and normal velocity. The variables F, g, Z, and C were de- 
fined as: 


F = u/U, 


g = H/Hg 


VS-,e 


(18) 

(19) 

(20) 


and 


C = py/PgUg 

Boundary conditions for the above system of equations were: 
at n = 0 

F = 0 


( 21 ) 


g = H^Hg (non-adiabatic) 
ag/an = 0 (adiabatic) 


V = V 


w 


6 



(porous wall) 


^12 

7 = 2 = —L£ 

w 


3y 


W 


3 Z/3n = 0 (slot injection) 


and at n * Hg 


F = 1 

g = 1 

Z = 1 

2.2 Parabolic Form of the Conservation Equations 

Equations 13, 14, 15, and 16 when expressed in the general pa- 
rabolic form become: 

Momentum ; 

F" + A^F‘ + A^F + A 3 + A^F^ = 0 



A2 ~ “BF/Aq 
A 3 = ePg/pAQ 

A 4 = -25F/Aq 
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and 


Aq = C(1 + e"") 


^0 " ^ 0 ^^ 

Energy : 

g" + A^g' + A2g + A3 + A^g^ = 0 

where 



A 


2 


= 0 




^ Pr^ ■ ^0 






and 


Species : 


where 


A4 = -25F/Aq 


A = — (1 + -•**— e^) 
^0 Pr Pr* ^ ^ 


Aq Aq/C 


A' + A^Z' + A2Z + A3 + A^Z^ = 0 


A = — 

ni — r- 


A ' 
C ' ^^0 


^ l^n 


A, 


A 2 - 0 
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^0 “ Pr Pr^. ^ ^ 

and 

^0 = 

2.3 Numerical Solution Procedure 

The finite-difference scheme used to solve the transformed equa- 
tions was an implicit method of the Crank-Nicolson type. The method has 
been used sucessfully by a number of authors including Anderson and Lewis 
(Ref. 11), Miner, Anderson, and Lewis (Ref. 12), Davis (Ref. 13), Blottner 
(Ref. 14), and Harris (Ref. 15). The equations were written in the gen- 
eral parabolic form 

W" + A^W + AgW + A3 + A^W^ = 0 (22) 

where W was the dependent variable F, g, or Z and the coefficients were 
functions of 5, n, W and W. The boundary layer was considered as a grid 
of nodal points with a varying step size in the normal direction as shown 
below. 
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Values of the dependent variables at m + 1 were calculated across 
the boundary layer beginning with the boundary conditions at n = n_ and 
using the relationship (Ref. 11) 




+ F. 


where 


-C 


^n " B 


„ + A E , 
n n n- 1 


F 


n 


Bn + K 

n n 


F 

E 


nd 

n-1 


The coefficients A, B, C, and D were evaluated from the coefficients of 
Equation 22 and the step sizes at each grid point using the following re- 
lations: 


A 


n 


■ ^In 

(An^ + 


(23) 


and 


. -2 ^ (''"n ' '^Vl* . . . ^ (24) 

''"n-1 


* ''in "Vl 
T% + in 

(25) 

. ^'In ''m,n 

(26) 

n AC 



The values of E-j and F^ were determined by the boundary condition at 
n = 0. For example, when the boundary condition set 3W/9n equal to zero, 
E^ became unity and F-j became zero. In another case where W at the wall 
was a given value, as with porous injection or with a non-adiabatic wall, 
then E.| was zero and F^ was the value of W at the wall. 
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When the equations are written in the general parabolic form, they 
form a coupled, non-linear system which must be solved by iteration. At 
each value of 5 , the species, energy, and momentum conservation equations 
were iterated using the above procedure and the continuity equation was 
solved by integration of the expression 

n 

V = J (2CF^ + F) dc (27) 

Iteration continued until Z, g, and F at each grid point across the bound- 
ary layer changed between iterations by less than the desired amount. 


2.4 Fluid Properties 

The fluid properties for air injected into air were obtained from 
the following standard relations. Density was computed from the state 
equation 

P = pM^/RT (28) 

and the viscosity was computed by Sutherland's formula 


y 

*^ref 


T X + C’ 
ref 


T + C* 



(29) 


where C* = 110.3°K. 

The transport properties were defined by setting Pr = 0.71, Pr^ = 0.9, 

Le = 1.0, and Le^. = 1 .0. 

For the case of foreign gas injection such as helium into air, ther- 
modynamic and transport properties were supplied by the tabulated data of 
Jaffe, Lind, and Smith (Ref. 16). Table 1 gives the data used to compute 
the properties of the individual species. The species specific heat at 
constant pressure was Cp^. = 31025 ft^/sec^-°R and then enthalpy was calcu- 
lated by h^ = Cp^ T ft^/sec^. The species specific heat at constant vol- 
ume was calculated from: 
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The mixture density was obtained as follows; 


where 



-Jlu— 


slug/ft^ 


Mi = = 4.0026 


Mf = M^^^ = 28.966 

The enthalpy and specific heat of the mixture were calculated by the 
following expressions: 


h = (l-C^)h^ + Cfh^ 


Cp = (l-C^)Cp^ + C^Cpf 

The mixture viscosity was calculated by Wilke's formula: 


^ = 1 + 


1 ^ G,.X./X^ 


(30) 


where 


and 


G 


Va 


X. = (1-C^)/M. [C^/M^ + d-C^)/M.] 
Xf = Cf/M^ [C^/M^ + d-C^)/M.] 


1 + 


M. 


\-l/2 



/^v\ 

rw 



2 


The mixture thermal conductivity was also calculated using Wilke's formula 
(30) where the species conductivity was calculated by the following expression 
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2.5 Eddy Viscosity Model 


The eddy viscosity, e , was computed by either a two-layer model 
consisting of an inner and outer region or a multi-layer model proposed 
by Beckwith and Bushnell (Ref. 2). Both models were based on Prandtl's 
mixing length concept given by 


* * *2 

e = p A 


3U 

ay* 


(31) 


The two-layer model used the eddy viscosity law of Van Driest for the 
inner law and an outer law based on Klebanoff's modification of Clauser's 

ic 

law (Ref. 17) while the multi-layer model varied the definition of £ in 
each layer. 


2.5.1 Two-Layer Model 

The two-layer model separates the boundary layer into two regions, 
one near the wall and the other consisting of the remainder of the bound- 

ic 

ary layer. A different expression for a in equation 31 is used in the 
two regions. In the present work, the expression developed by Van Driest 

(Ref. 18) was used in the wall region 

/ = k^ y* [l - exp (-yW)] (32) 

where k-j = 0.4 

^ * * * 
y = y /v 




(t„ /p ) 


1/2 


For a non-porous wall without a pressure gradient, = 26. For 
a case with mass transfer through the wall, A^ was corrected as proposed 
by Cebeci (Ref. 19) 


A"*" = 26 exp (-5.9 


(33) 
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where 


''w ' \ /“f 

For cases with both mass transfer and pressure gradients A"*” becomes 
(Ref. n) 

A* = [exp (11.8 V ■'■) - l] + exp (11.8 v '^) \ 

I I- W J W j (3,) 


^ v*/p* (u/)3 
dx ^ 


The eddy viscosity law in the region near the wall becomes 
" 4 = (k, y*)^ [1 - exp (V/A-^)] 2 < 4 

p ' L J 3y p 

For the outer regions of the boundary layer, the eddy viscosity was 
based on the work of Clauser and is given by 

+* * * ★ * 
e = kg p Ue 

* 

where kg = 0.0168, 6|^ is the incompressible, two-dimensional boundary- 
layer displacement thickness and y is the Klebanoff intermittency factor 
which was approximated by 

Y = [l + 5.5 (y*/6*)® ] (37) 

In the two-layer model the outer law was first used across the entire 
boundary layer. Then, starting at the wall, the inner law was used up througf 
the boundary layer until the eddy viscosity from the inner law was greater 
than or equal to that of the outer law. 


2.5.2 Multi-Layer Model 

The multi-layer eddy viscosity model is due to Beckwith and Bushnell 
(Ref. 2) and also used by Miner and Lewis (Ref. 4), The Prandtl mixing 
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length concept was modified using the Van Driest damping factor to 
give the following expression for in dimensional variables 


[l - exp (V/A'^)] 

The mixing length £ Is then calculated in the different regions as follows 

For a fully-developed turbulent boundary layer, the flow is divided 
into three layers by specific values of y and the mixing length is de- 
fined in each of these layers as: 


9u (38) 

sy y 


Point 

0 

1 

2 


y/6 

0 

0.1 

0.3 


where K = 0.4 and f is given by 


0 

0.1 K6 
f6 


f = 0.265 

★ ^ 

where /6|^ and 6|^ is 

thickness. 


- 0.196 Hj,* + 0.0438 

the incompressible boundary-layer momentum 


For tangential slot injection, the above model was modified to con- 
sider three distinct regions or zones (Fig. 1). These zones were deter- 
mined by the relative values of the mixing lengths in the slot flow re- 
gion, in the mixing region, and in the outer region of the bound- 

ary layer, 

For the slot flow region the mixing length was constant and given by 


where s is the slot height and a^ is a constant. The suggested value for 

turbulent slot flow was a =0.14 (Ref. 2). 

s 

In the mixing region, was given by 

“ 3m W Pr./Le+ 
mm t t 
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where W was the width of the mixing zone between the slot flow and 
the outer flow and was a constant. The values of should range from 
0.05 to 0.12, with 0.09 being the recommended value (Ref. 2). The 
mixing zone width W is defined by 


w = yf - y„ 


where y^ = y at Z * 0.99, = y Z = 0.01 and Z is given by 


= (C, 


c. )/(c. - C. ) 

1 ' 1 ,e 1 ,w^ 


The mixing length in the outer flow was defined as 


‘b ' ^ - y„> 

where f^ was 

f^ = 0.265 - 0.196 H* ^ + 0.0438 H* ^ 

S K 5 S K 5 S 

and Hi, _ = Sk ./e^ _. The subscript s denotes that the lower limit of inte- 

K)S KjS KjS 
it 

gration for and was changed from y = 0 to y = as recommended by 
Beckwith and Bushnell (Ref. 2). 

With the above definitions of and £^, the three zones may 

be specified. Zone 1 is the initial mixing zone and is defined by the 

inequal ity 




m b 


The model used in Zone 1 was a five-layer model using the following 
pivot points and coordinates of y and z: 
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Point 

X. 

_£ 

0 

0 

0 

1 

wh 

1 ®s 

2 


1 ®s 

3 

^m 

^ a W 
Le^ m 

4 

‘ * r ) 1 



where s is the slot height, t is the lip thickness, = y at Z = 0.5, 

is the thickness of the boundary layer above the lip, and 

6,^ = 6^ + s + t. 
r 0 

Zone 2 was the intermediate mixing zone and used a four-layer model 
formed by dropping pivot point 2. This zone was defined by the inequality 

s — m — b 


and the pivot points and coordinates of y and a are 


Point 

0 

1 


0 


Le^K 


0 

3m w 
Le^ m 


3 


4 


'm 


rr- 3m w 
Le m 




% 
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Zone 3 approached the equilibrium boundary-layer state and was given 


by 


> £ 

s — m — 


b 


This zone used the three layer model and the pivot points are given by 


Point 

0 

1 

4 


Z 

0 

0.1 


0.3 


£ 

0 

0.1 K6 


A schematic representation of the mixing length profiles, velocity pro- 
files, and species concentration profiles is presented in Figure 1. 


2.6 Application of the Finite-Difference Method 

In order to begin the numerical solution, initial profiles for the 
dependent variables Z, F, and g are needed. These guesses were of the 
following form: 


at n < ng 


at Hg 

F = 1 - e'^ 


F = 1 

F' = e"’^ 


' F' = 0 

F" = -e"'’ 


F" = 0 

g = Ml - 

9w^' 

9 = T 

g' = (1 - g„) 

F' 

g' =0 

g" Ml - g„) 

F" 

g" = 0 

Z = 1 


Z = 1 

Z' = 0 


Z' = 0 

r = 0 


Z" = 0 
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The solution of the continuity equation across the boundary layer 
was assumed as 


V = Vw - n 

and the temperature profile was calculated from the g profile 

U ^ U ^ « 

T/Tg = (1 + 1/2 g - 1/2 

e e 

u ^ u/ 

(T/Tg)' = (1 + 1/2 g' - j^FF' 

e e 

The density profile was assumed to be that of a perfect gas 


(39) 

(40) 



and the Chapman-Rubesin Factor, C = pw/pgPg was set to unity across the 
boundary layer. 


Starting with the above guesses, the equations were solved to deter- 
mine Z, g, and F. New values of T/T^, (T/T^)', p/p., V, and C were com- 

6 6 6 

puted where 


n 

V = - j (2cF^ + F) dr, 

0 


(41) 


C = 



where C = 


(1 

(t/Tg + C) 


(Sutherland's Law) 


or 

C (Power Law) 
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For tangential slot injection, the initial profiles for F, g, and 
Z must be supplied at the slot location. These may be obtained from ex- 
perimental data or by combining the results of a flat plate, finite- 
difference solution with an approximation for the slot profile. In the 
present work, flat-plate calculations were compared with experimental 

Mach number profile data to match as closely as possible the boundary- 

* 

layer thickness 6, the compressible displacement thickness 6 , the 
momentum thickness e, the Mach number profile and the velocity profile. 
The profiles which best matched these data were combined with an approx- 
imation for the slot flow and used as initial profiles at the slot 
location. 

Boundary-layer parameters were calculated from the converged so- 
lutions of g, F, and Z at each value of c* The definitions used to 
calculate these parameters were: 

Boundary-Layer Thickness: 

* 

The boundary- layer thickness, 6, was the value of y where 
u /Ug = 0.995 and was determined by interpolation of the F array. 


Incompressible Displacement Thickness : 


or 





1/2 


e e 




(1 - F) 


dn 


(42) 


(43) 
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Compressible Displacement Thickness : 


* 

6 = 



1 


* *"| 

P U 


* * 



* 

dy 


(44) 


or 


r K 

a* *^n 

e e 0 


dn 


(45) 


Momentum Thickness: 


* 

0 = 



* * 


P U 




(46) 


or 


* e '^e 

Vr -1 

a PflU r'^ 

^e e 


F (1-F) dr 


Skin Friction Coefficients: 


•Tw 


'f *2 

00 p U 

^ on 


where x^* = 


3U 

* 

3 y 


W 


or 


2e»in W« C Pa U_ r"^ 
VD e w e e 


(25) 


1/2 


3 £ 

3 n 


w 


(47) 


(48) 


(49) 
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2.7 Global Pressure Interaction 


The present finite-difference method contained, as an option, the 
iterative procedure developed by Miner and Lewis (Ref. 4) to include the 
effects of pressure interaction. The first boundary-layer calculation 
was made using the input pressure distribution dp/dx = 0 and calculated 
the displacement thickness distribution. The method then used the dis- 
placement thickness to calculate a new pressure distribution and began 
a new boundary-layer solution. This iteration was continued until the 
change in the pressure distribution between iterations was sufficiently 
small . 

The method assumed that the flow was isentropic and that Prandtl- 
Meyer theory was applicable. From the displacement thickness slope, 

d6*/dx, the edge Mach number was found using Prandtl -Meyer theory. Since 

* 

the 6 distribution may not be smooth, a six-point walking least squares 
log-log curve fit was used to calculate d6*/dx. This curve fit provides 
a smoother derivative than a three or four point Lagrangian interpolating 
polynomial. Isentropic relations were used to calculate p/p^ from the 
Mach number. Since the interaction between the pressure distribution and 
the displacement thickness may be rather strong, a weighting factor w was 
used so that 


Pnew = Pold Pi* 

where p was the pressure distribution to be used in the next iteration, 

^new ^ 

Poid was the distribution from the last iteration, and p^ was the dis- 
tribution computed from 6*. Suggested values of w range from 0.05 to 
0.2 depending on the degree of difficulty in maintaining stability. 
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III. RESULTS AND DISCUSSION 


The results of the present method were compared with a recent ex- 
perimental study conducted by Schetz and Van Overeem (Ref. 10). The 
experiments were conducted in the supersonic wind tunnel at Virginia 
Polytechnic Institute and State University. Air was injected through a 
0.635 cm (0.25 in.) slot and/or a 3.81 cm (1.5 in.) porous plate. The 
width of the wind tunnel test section was 22.86 cm (9.0 in.) which gave 
for the slot an injection area of 14.52 cm^ (0.01563 ft.^) and for the 
porous plate, an area of 82.26 cm (0.08854 ft. ). Since combinations 
of slot and porous geometries were considered, the results were present- 
ed in terms of the total mass flow m = p^u.A. and non-dimensional ized by 

J J J p 

a reference freestream value defined as p^u^ (1 ft. ). These experiments 
were conducted at = 2.9, with a stagnation temperature of 294°K, and 
a stagnation pressure of 6.803 atmospheres. Pitot pressure profiles, 
wall pressure, and wall shear stress data were taken throughout a range 
of mass injection rates. The nozzle was modified so that the streamline 
along the axis of the two-dimensional, symmetric nozzle was replaced by 
a solid surface (see Fig. 2a). Since the lower portion of the modified 
nozzle was a flat surface, a flat-plate geometry was used in the present 
method to simulate the experimental conditions. The use of this geometry 
caused some problems in determining the initial value of the stream- 
wise coordinate. The nozzle boundary layer begins forming somewhere up- 
stream of the curved portion of the converging section of the nozzle 
while the flat-plate boundary layer is assumed to begin at the leading 
edge of a constant pressure flat plate. Thus the virtual flat-plate 
length had to be determined to best simulate the boundary layer on the 
nozzle wall. Experimental pitot pressure profiles, taken at the slot 
location, were compared with calculations of the flow over a flat plate. 
The predicted profiles best matched the measured profiles at a distance 
of 53.34 cm (21.0 in.) downstream of the leading edge. Figure 3 shows 
good agreement between the flat-plate prediction and the experiment at 
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this location. The value of 5 = 7.8975 at 53.34 cm was used as the 
initial value of c for the mass injection study. 

Initial profiles were taken from the experimental data of Ref. 10 
at the slot location. Pitot profiles were measured at the slot and used 
to obtain Mach number profiles, assuming a constant static pressure. By 
assuming a constant total enthalpy profile (H/H^ = 1.0), initial velocity 
profiles were calculated from the Mach numbers. The initial air species 
profile for helium injection was set to unity above the slot lip and to 
zero below the lip. For the cases using air injection the species pro- 
file was set to unity across the boundary layer. 

Flat-Plate Results 

The above determined initial profiles were used to simulate a boundary- 
layer flow over a flat plate. The predicted skin friction was 40% higher 
than measured by the floating-element balance. No Preston tube measure- 
ments were obtained on the solid flat-plate; however, extrapolation of the 
Preston data to m = 0 presented in Fig. 9, gave considerably higher values 
of wall shear than were measured by the floating-element balance. To help 
resolve these differences, the results were compared with previous investi- 
gations. Figure 4 compares the skin friction under the present conditions 
with data by Coles (Ref. 20) and the numerical results of Anderson and Lewis 
(Ref. 11). The present results are in excellent agreement with Anderson 
and Lewis and fall between Cole's Case 26 and Case 20. The extrapolated 
Preston tube value also agrees well with the present method but the balance 
data on both the solid plate and the porous plate with m = 0 were low. 

Since the skin-friction distribution was not measured with either the 
balance or the Preston tube, it was impossible to resolve the differences 
between the balance data and the present prediction. Figure 5 shows the 
measured velocity profile at the slot location in law-of-the-wall coor- 
dinates. Again the numerical results and extrapolated Preston tube data 
are in good agreement while the balance data are below the prediction. 
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without skin-friction distributions it was impossible to determine whether 
or not the balance and Preston tube data show the same flat-plate trends. 

Figures 6-8 show the Mach number profiles at three locations down- 
stream of the slot. 

Normal Injection Through a Porous Wall 

The first mass transfer geometry considered was a porous wall. 
Porous-wall injection caused a large reduction in the skin friction, but 
the effectiveness of the porous wall injection was limited by boundary- 
layer separation. Figure 9 presents the skin-friction variation as a 
function of the mass transfer rate. The balance and Preston tube data 
again do not agree. At high injection rates, the balance data were 
slightly lower than the Preston tube data but had a similar slope. At 
low injection rates, the agreement was poorer. The limited Preston tube 
data indicated a much more rapid reduction in skin friction with mass 
transfer than the balance data, and the Preston tube predicted a higher 
skin-friction at m = 0. Differences in the two measurement techniques 
create some question as to the effectiveness of normal injection. The 
Preston tube data show the porous wall was the most effective geometry per 
unit mass in reducing skin friction. The balance data show it was the 
least effective. More experimental data are needed to determine conclusively 
the effectiveness of the porous wall. The numerical study showed very 
good agreement with the Preston tube data and predicted boundary-layer 
separation over the porous plate when the mass-transfer rate was greater 
than m/p^u^ (1) = 0.0016 or Py^v^p^u^ = 0.01807. This separation would 
account for the sudden loss in effectiveness evident in the Preston tube 
data. The balance data show the porous-wall technique suffered from a 
roughness-induced increase in skin friction, but no Preston tube data 
were available to confirm this increase. Since the balance did measure 
lower skin-friction values on the flat plate, this apparent increase 
should be confirmed by additional experimental data. It is probable 
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that the balance was more sensitive to roughness effects since it was 
flush with the wall and completely enveloped by the sublayer. The 
Preston tube, on the other hand, could not be placed on the surface due 
to the finite thickness of the tube walls and may be above the sublayer 
which was most affected by the roughness. More Preston tube and balance 
data are needed to better determine the effects of roughness of the porous 
surface on the skin friction. 

Calculations were also made which included the effects of pressure 
interaction. The pressure distribution resulted in higher skin-friction 
values than did the case with a constant edge pressure. As the mass in- 
jection rate increased, the effects of pressure interaction increased 
because the predicted pressure gradients became stronger over the porous 
section (see Fig. 10). These increased pressure gradients also caused 
the boundary layer to separate at lower mass injection rates than the 
cases without pressure interaction. Mach number profile comparisons at 
two different Injection rates are presented in Figs. 11-16. At the lower 
injection rate, the agreement with experimental measurements was excellent 
and pressure interaction had very little effect. At the higher injec- 
tion rate, the agreement was rather poor with the predicted Mach numbers 
being as much as 20% lower than the experimental data. Pressure inter- 
action solutions showed better agreement but were still as much as 10% too 
low. As the external stream flowed over the porous plate, normal injec- 
tion caused a thickening of the boundary layer. Downstream of the porous 
section, the boundary layer became thinner after the normal injection 
stopped. In the low injection case, the thickening of the boundary layer 
was not very severe and the flow recovered quickly thus giving good agree- 
ment with the profiles taken downstream of the porous section. At the 
high injection rate, however, the boundary layer became much thicker and 
could not adjust as quickly; thus the present method predicted a much 
thicker boundary layer downstream of the porous wall. The pressure gra- 
dient behind the porous section was favorable and resulted in better 
agreement with the experimental profiles. 
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Pitot pressure profiles for the first two stations at both injec- 
tion rates are presented in Figs. 17-20. The results were similar to 
those of the Mach number profiles except near the outer edge of the bound 
ary layer at station one. Here the experimental data do not asymptot- 
ically approach unity at the edge. This condition was caused by shock 
waves generated at the junction of the injection box and the test sec- 
tion wall as shown by the Schlieren photographs in Ref. 10. At station 
two, the probes were completely behind these waves and the profiles did 
not show a jump near the outer edge of the boundary layer. 

Tangential Injection Through a Slot 

The second injection geometry consicfered was a tangential slot. 

In this part of the study, both eddy viscosity models were used and a 
foreign gas injectant was briefly considered. Figure 21 presents the 
skin-friction variation with the mass-transfer rate. The balance data 
were again lower than the Preston tube data, but the experimental re- 
sults did show the same trends. The Van Driest-Clauser eddy viscosity 
model predicted a skin-friction distribution about 30% higher than did 
the Beckwith-Bushnell model. The difference is due to the differences 
in the velocity profiles caused by changing eddy viscosity models (see 
Fig. 28). The Beckwith-Bushnell model agreed well with the Preston 
tube data. Pressure interaction did not have the effect on slot injec- 
tion that it had on the porous wall geometry since the pressure gra- 
dients were weaker as shown in Fig. 22. The presence of the step in 
the wall caused a 12% reduction in the skin friction from the flat- 
plate value as shown in Fig. 4. The skin friction continued to decrease 
with an increase in mass transfer until the matched pressure injection 
condition was reached near m/p„u^ (1) = 0.0016. This mass flow rate 
corresponds to X = p,-U 4 /p^u^ = 0.1024, and the slot Mach number was 
approximately 0.6. As the mass flow rate continued to increase, the 
turbulence level in the slot became high enough to cause an increase in 
the wall shear stress. 
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Calculations were made to compare the wall pressure distribution 
with the experimental data. For this part of the study, the numerical 
method used the same edge pressure as the experiment, Pg/Pg = 0.031. 

This resulted in a freestream Mach number of 2.91 which was slightly 
higher than the 2.8 used in the skin-friction study. Figure 22 shows 
the predicted wall pressure distribution was in reasonable agreement 
with the measured values. 

A special case using helium as the injectant was briefly considered. 
Figure 23 shows the skin-friction variation with mass transfer using the 
two-layer eddy viscosity model. The agreement with the experimental data 
was good. Comparison of Figs. 21 and 23 shows that helium was more effec- 
tive in reducing the skin friction than air. No other experimental data 
are available from the helium injection study. 

Comparisons of the predicted and experimental Mach number profiles 
for the slot configuration are made in Figs. 24-26 for the case of 
m/p^u^ (1) = 0.00072. Figures 24-26 compare the profiles from both 
eddy viscosity models with the experimental results at m/p^u^ (1) = 0.00072. 
At station 1, the Van Driest-Clauser model gave excellent agreement in 
the outer region of the boundary layer but did not agree well near the 
wall. The Beckwith-Bushnell model gave better agreement throughout the 
entire boundary layer especially near the wall. Pressure interaction 
had little effect on the multi-layer model but resulted in much better 
agreement when the two-layer model was used. Figure 22 shows the pres- 
sure was recompressed to the freestream value more quickly using the 
Van Driest model than when using multi-layer model. The pressure gra- 
dient at the first probe station was larger using the Van Driest model 
and thus had a greater effect on the solution. At station 2, the 
Beckwith-Bushnell model again agreed better near the wall than did the 
two-layer model. Here pressure interaction had a greater effect on the 
multi-layer model because the Van Driest model had recovered to the free- 
stream pressure and experienced a small pressure gradient. A weak adverse 
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pressure gradient existed downstream of the slot as the flow first over- 
expanded around the step and then was recompressed to the freestream 
pressure causing a thicker boundary layer to be predicted by the pres- 
sure interaction solution. Profiles at the third station showed the 
two-layer model results were in better agreement with the data than the 
multi-layer model. The pitot profiles presented in Figs. 27-29 show 
the same agreement as the Mach number profiles. 

In general, the Beckwith-Bushnell eddy viscosity model gave better 
agreement near the slot, while further downstream, the Van Driest-Clauser 
model more closely matched the experimental data. The present method 
underpredicted the Mach number profile considerably. This was due to 
the low mass injection rate which corresponds to a = 0.046 and which 
was about 40% lower than the low pressure injection rate used in Refs. 

5 and 6. Miner and Lewis (Ref. 6) found the present method overpredicted 
the Mach number profiles at matched pressure injection but gave good 
agreement at the low pressure rate. The present work shows that at even 
lower injection rates the method underpredicts the Mach number profiles. 

Combined Slot and Porous Wall Injection 

A combination of the two previous geometries was considered in the 
third part of this study. Air was injected through the slot and porous 
plate at varying mass-flow rates. Calculations predicted separation 
would occur over the porous section when more than 30% of the total mass 
flow was injected through the porous wall. Figure 30 compares the skin- 
friction reduction using both eddy viscosity models for the cases of 100% 
slot injection and 75% slot/25% porous injection. The Beckwith-Bushnell 
model predicted lower skin-friction than the Van Driest-Clauser model 
because of differences in the velocity profiles near the slot as seen 
in Figure 32. 

One would expect combined injection to show characteristics of both 
the porous wall and slot injection geometries. Normal injection through 
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a porous wall decreased the skin friction markedly but separation occurred 
at mass flow rates of m/p^^u^ (1) = 0.0017. Slot injection also reduced 
the skin friction but reached a mass flow rate where the skin friction 
reached a minimum and then began to increase. The calculations for com- 
bined injection did indeed show characteristics of both geometries. At 
low mass injection rates, the porous plate reduced even further the skin 
friction downstream of the slot, but at the higher injection rates, the 
slot turbulence level increased causing the skin friction to increase. 

The most important effect of downstream normal injection was the decrease 
in the minimum skin friction attainable. The Beckwith-Bushnell model 
predicted a 50^ lower skin friction with the 75% slot/25% porous geometry 
than with the slot-only case. • The minimum was shifted to higher injection 
rates for combined injection because only a portion of the total mass flow 
was transferred through the slot, and thus the slot Mach number was lower 
even though the total mass flow was the same. 

Pressure interaction had little effect on the slot flow but in- 
creased considerably the skin friction on the porous wall. A similar 
result was evident in the combined slot/porous injection case. At low 
injection rates, where the normal injection was most effective, pressure 
interaction solutions predicted higher skin-friction values. At higher 
rates, where the slot injection effects were stronger, pressure inter- 
action caused only small differences. 

Insufficient experimental data were available to validate the pres- 
ent predictions and most of those were from the skin-friction balance. 

The presence of the porous section again caused a large increase in the 
shear stress measured by the balance. At a mass flow of m/p u (1) = 0.00074, 
the balance measured a 76% increase over the slot-only configuration. 

As mentioned above, the balance was probably more sensitive to surface 
roughness than was the Preston tube, but it remains to be determined if 
roughness alone could account for such a large change. The Preston tube 
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measurements were not complete enough to determine whether or not sur- 
face roughness affected the Preston tube. The differences between the 
prediction and experiment could not be resolved with the limited ex- 
perimental data available. It is impossible to completely validate the 
predictions from the numerical method until more experimental data are 
available. 

Predictions of the wall pressure distribution again showed only rea- 
sonable agreement with the experimental data as shown in Fig. 31. Mach 
number profiles shown in Figs. 32-34 and pitot pressure profiles shown 
in Figs. 35-37 were similar to results for the slot geometry. The 
Beckwith-Bushnell eddy viscosity model again showed better agreement 
near the slot than did the Van Driest-Clauser eddy viscosity model. Fur- 
ther downstream, however, the two-layer model more closely matched the 
experimental data than did the multi-layer model, but both eddy viscos- 
ity laws predicted a thicker boundary layer than was measured. This 
was caused by the high percentage of the flow entering through the por- 
ous wall. When more than 30% of the flow was injected through the wall, 
separation was predicted to occur over the porous section. The 70% 
slot/30% porous injection case was very close to this limit, and there- 
fore the flow was near blow-off over the plate. The flow did not re- 
cover quickly downstream of the porous wall, and thicker profiles were 
predicted at the measurement stations. 

Pressure Interaction solutions using the two different eddy viscos- 
ity models showed completely reversed trends. Using the Van Driest model, 
pressure interaction caused a thicker boundary layer, but using the 
Beckwith-Bushnell model, a thinner boundary layer was predicted. Figure 
31 shows the edge pressure was actually decreasing behind the porous 
plate using the multi-layer model. This favorable pressure gradient 
caused the thinner boundary layer while the solution using the Van Driest 
model was subjected to an adverse pressure gradient which thickened the 
boundary layer. 
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Figure 38 shows a comparison of the skin-friction variation with the 
mass flow rate for the three geometries considered. According to the pres- 
ent predictions, the porous wall geometry provided the greatest reduction 
in skin friction per unit mass of injectant but was limited to low injec- 
tion rates by boundary-layer separation. Tangential slot and combined in- 
jection geometries both produced a reduction of wall shear stress due to 
the presence of the step in the wall. Slot injection decreased the skin 
friction as the mass flow increased until the turbulence level became 
great enough to cause an increase in the wall shear stress. The combined 
injection model predicted even lower skin friction at low flow rates than 
did the slot geometry. The greatest effect of the combined injection was 
to reduce the minimum in skin friction at the matched pressure slot in- 
jection condition caused by injecting a portion of the flow through the 
wall. This delay allowed the slot/porous injection to produce a lower 
minimum skin friction before the slot turbulence became large enough to 
increase the wall shear stress. The numerical results and Preston tube 
data show normal injection was the most effective in reducing skin friction 
per unit mass injected. Schetz and Van Overeem (Ref. 10) concluded that 
the slot injection was the most effective based on the balance data. The 
reason for this contradiction was the balance data indicated normal injec- 
tion to be far less effective in reducing skin friction than the Preston 
tube data or numerical results predicted. More data are needed for compari- 
son to resolve the differences between the Preston tube, floating-element 
balance, and the present predictions. 
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IV. CONCLUSIONS 


In the present study a finite-difference method was used to predict 
skin-friction reduction for porous wall and/or slot injection geome- 
tries. A comparison of the three injection geometries showed that nor- 
mal injection provided the greatest reduction in skin friction per unit 
mass of injectant but was limited by boundary-layer separation. The 
tangential slot and combined injection geometries both produced a lower 
skin friction at m * 0 due to the step in the wall. As the mass injec- 
tion through the slot increased, the skin friction decreased until the 
mass flow rate became so large that the turbulence level caused the skin 
friction to increase. 

Combined injection offered some improvement over the slot geometry 
at low mass injection and showed considerable promise at the higher injec- 
tion rates. With a portion of the mass flow coming through the porous 
wall, the increase in slot turbulence and skin friction was delayed. 

This situation resulted in predicting a much lower skin friction before 
the slot turbulence became great enough to increase the skin friction. 

The numerical calculations predicted separation would occur over the 
porous section of the wall when more than 30% of the total mass flow was 
injected through the porous wall. Comparisons of the predictions with 
the limited Preston tube data were good while comparisons with the balance 
data were poorer. Differences exist in the data from the Preston tube 
and the floating-element balance. First of all, the balance data were 
lower than previous zero pressure-gradient flat-plate investigations at 
similar Mach numbers while no Preston tube measurements were available 
from Schetz and Van Overeem. The balance also measured an increase in 
skin friction over the flat plate when the porous section was added. 

This increase was attributed to surface roughness by Schetz and Van Overeem 
(Ref. 10) while no Preston tube data were taken for comparison. The 
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balance data indicate that normal injection was not the most effective 
.injection geometry, whereas the limited Preston tube data show just the 
opposite. For the combined slot/porous geometry, the balance measured 
a large increase in skin friction over the slot-only case, but the 
Preston tube data showed no such increase. More experimental data are 
needed to resolve the differences in the two measurement techniques and 
to determine the effects of surface roughness on the skin friction. 

Mach number profile comparisons showed good to excellent agreement 
with the experimental data at low injection rates. At the higher mass- 
transfer rates, the present calculations predicted a boundary-layer 
thickness that was almost 25 % greater than the measured value. For the 
slot and combined geometries, the Beckwith-Bushnell eddy viscosity model 
gave better agreement near the slot than did the Van Driest-Clauser 
model. Further downstream, the two-layer model more closely matched the 
experimental data than did the multi-layer model. 

Pressure interaction caused a substantial increase in the porous 
plate skin friction but made little difference in the slot injection 
case. With combined injection, pressure interaction increased wall shear 
at low mass flow rates but had little effect at the higher rates. The 
global pressure interaction scheme also gave a reasonable prediction of 
the wall pressure distribution. On the porous plate, the favorable 
pressure gradient caused a thinner boundary layer in the pressure inter- 
action solutions. For the slot injection a weak adverse pressure gradient 
resulted in a thicker boundary layer. 
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Table 1 

g ■ ■■ ~ — — 


Thermodynamic and Transport Property Data for Helium Injection 


Air Specific Heat at Constant Pressure 


Cp = A + BT + CT^ + DT^ + ET^ + FT^ 


ft^/sec^ - °F 


A = 6.03517 (10’^) 

B = -9.45091 (10"^) 


0 R < T £ 2000 

C = -7.30226 (10'^) 
D = 1.73227 do"®) 


E = -9.76574 (10’^°) 
F = 1.74651 (10‘^^) 


Air Static Enthalpy 
hj = AT + BT^ + CT^ + DT^ + ET® + FT® 
0 °R < T < 2000 °R 

A = 6.03517 (10^) C = -2.43408 (10"^) 

B = -4.72545 (10"^) D = 4.33069 (10“^) 


2 2 
ftVsec^ 


E = -1.95314 (10'^°) 
F = 2.91086 (10"^^) 


yj X 10^ = A + BT 


A = -1.93368 (10"^) 
B = 8.93259 (10"^) 

A = 5.50269 (10"^) 

B = 7.94717 (10"^) 


Viscosity 

CT^ + DT^ + ET^ + FT® 
90 °R < T < 6300 °R 

Air 

C = -3.9185 (10“®) 

D = 1.17285 (10'^) 

Hel ium 

C = -2.21781 (10“®) 

D = 4.04477 (10'^°) 


Ib^ - sec/ft^ 


E = -1.64934 (10"^^) 
F = 8.75154 (10'^®) 

E = -1.80635 (10’^'^) 
F = -1.01127 (10'^®) 


Binary Diffusion Coefficient 
p D^.j = A + BT + CT^ + DT^ + ET^ + FT® Ib^/sec 
90 °R < T < 6300 °R 


A = -9.75983 (10"^) 
B = 1.43997 (10‘®) 


Hel ium 

C = 3.99118 (10"®) 

D = -4.05169 (10‘^°) 


E = 4.70538 (10"^^) 

F = -2.35942 (10"^®) 
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VELOCITY PROFILES 



MIXING LENGTH PROFILES 



SPECIES CONCENTRATION PROFILES 

Figure 1: Schematic of Slot Injection Flow Field 
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Slot Location 




b) NUMERICAL STUDY 

Figure 2: Schematic of Experimental Nozzle and Model for 

Numerical Calculations 
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Experiment of Schetz and Van Overeem 
Porous Flat Plate with No Injection 

Extrapolated Preston Tube Data ID • 0. 241 cm 
Floating Element Balance ▲ 


Solid Flat Plate 

Floating Element Balance # 



Figure 4: Flat Plate Skin Friction 




Experiment of Schetz and Van Overeem 

• Solid Flat Plate-Floating Element ■ Porous Flat Plate-Preston Tube D • 0. 241 cm 

APorous Flat Plate-Floating Element m-0 

O Cole's Case 20 IVl^o -3.701 Rex-3.54xloJ 

Van Driest (1956) -3.701 Rey -3.54 x lOj 

P resent Method Mo, - 2. 8 Re^ - 3. 86 x 10 ' 



Figure 5; Velocity Profiles at Slot Location, Law of the Wall 
Variables 
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Results of Present Method 
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Figure 6: Flat Plate Mach Profiles, Station 1 
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Results of Present Method 
— Van Driest -Clauser Model 
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•1 Probe-^ 
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294 °K(530 °R) 


6.803 atm (100 psia) 

293.8 '"K(529 °R! 
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Mach Number 


Figure 8: Flat Plate Mach Profiles, Station 3 






X 10 


POROUS INJECTION 


Experiment of Schetz and Van Overeem 
• Balance 

APreston tube, D « 0.241 cm 
Present Method 
Van Driest-Clauser 

Van Driest-Clauser with pressure interaction 


M„ • 2. 8 Tq - 294 ® K (530 ® R) Pq - 6. 803 atm ( 100 psia) 



0.00 0.01 0.02 0.03 0.04 

p V /p u 
W W * “ 


Figure 9: Skin Friction Coefficient Variation with Mass Flow Rate, 

Porous Injection, Station 1 
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POROUS WALL 
Results of Present Method 


m/PooUoo (1) -5.52x 10"^, Van Driest 


m/p«,u„ (1) - 10.34 X 10"^, Van Driest £+ 
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X (cm) 

Figure 10: Predicted Wall Pressure Distribution for Porous Plate 

Injection 
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POROUS PLATE 
m/pa,Uoo (1) “5.5 X 10 ^ 
Station 1 

• Experiment of Schetz and Van Overeem 
Results of Present Method 
Van Driest-Clauser e*'’ Model 


Van Driest-Clauser with pressure interaction 



Mach Number 

Figure 11: Porous Injection Mach Profiles, Low Mass Injection, 

Station 1 
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Figure 12: Porous Injection Mach Profiles, Low Mass Injection, 

Station 2 
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Figure 13: Porous Injection Mach Profiles, Low Mass Injection, 

Station 3 
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POROUS PLATE 
m/p„u„ (D-10.4 X 10“^ 
Station 1 

• Experiment of Schetz and Van Overeem 
Results of Present Method 
Van Driest-Clauser Model 


- Van Driest-Clauser with pressure Interaction 



Mach Number 


Figure 14: Porous Injection Mach Profiles, High Mass Injection, 

Station 1 
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Figure 15: Porous Injection Mach Profiles, High Mass Injection, 

Station 2 
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Figure 16: Porous Injection Mach Profiles, High Mass Injection, 

Station 3 
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Figure 17: Porous Injection Pitot Pressure Profiles, Low Mass 

Injection, Station 1 
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Figure 19: Porous Injection Pitot Pressure Profiles, High Mass 

Injection, Station 1 
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Figure 20: Porous Injection Pitot Pressure Profiles, High Mass 

Injection, Station 2 
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Figure 21: Skin Friction Coefficient Variation with Mass Flow 

Rate, Slot Injection, Station 2 
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SLOT INJECTION 


• Experiment of Schetz and Van Overeem 
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Results of Present Method 
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Figure 22 : Wall Pressure Distribution, Slot Injection 
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SLOT INJECTION-HELIUM 


• Experiment of Schetz and Van Overeem-Balance Data 
Present Method -Van Driest 
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Figure 23: Skin Friction Coefficient Variation with Mass Flow Rate, 

Helium Slot Injection, Station 2 
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Figure 24: Slot Injection Mach Profiles, Low Mass Injection 

Station 1 
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Figure 25: Slot Injection Mach Profiles, Low Mass Injection. 

Station 2 


62 


Y (cm) 




Y (inches) 


SLOT INJECTION 

m/p«,Uoo (1) -7.2x 10"^ 0.046 

Station 3 

• Experiment of Schetz and Van Overeem 
Results of Present Method 
£■*■ No Pressure Effects Pressure Effects 


Van Driest-Clauser 

Beckwith-Bushnell 



Figure 26: Slot Injection Mach Profiles, Low Mass Injection, 

Station 3 
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Figure 27: Slot Injection Pitot Pressure Profiles, Low Mass 

Injection, Station 1 
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Figure 28: Slot Injection Pitot Pressure Profiles, Low Mass 

Injection, Station 2 
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Figure 29 : Slot Injection Pitot Pressure Profiles, Low Mass 

Injection, Station 3 
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Figure 30: Skin Friction Coefficient Variation with Mass Flow 

Rate, Combined Injection, Station 2 
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Figure 31: Wall Pressure Distribution, Combined Injection 
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Figure 33: Combined Injection Mach Profiles, Station 2 
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Figure 34: Combined Injection Mach Profiles, Station 3 
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Figure 35: Combined Injection Pitot Pressure Profiles, Station 1 


72 


Y (cm) 




Y (inches) 


COMBINED INJECTION 

ni/PcoUeod) “ 1.324 X 10"^ 70% Slot/30% Porous 

Station 2 

•Experiment of Schetz and Van Overeem 
Results of Present Method 

Van Driest-Clauser £■*■ 

Beckwith-Bushnell e+ 



Figure 36: Combined Injection Pitot Pressure Profiles, Station 2 
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Figure 37: Combined Injection Pitot Pressure Profiles, Station 3 
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Figure 38: Comparison of Porous, Slot, and Combined Injection Skin 

Friction Variation with Mass Flow Rate 
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